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Abstract 

I give a brief overview of the status of lattice QCD, concentrating on topics 
relevant to phenomenology. I discuss the calculation of the light quark spectrum, 
the lattice prediction of ay^{Mz), and the calculation of fs- 
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1 INTRODUCTION 



Lattice gauge theory has been somewhat out of the mainstream of particle physics for the 
past decade. It seems to me, however, that the field is now coming of age. It has certainly 
grown rapidly: roughly 300 people attended the recent LATTICE '92 conference, compared 
to about 130 at LATTICE '86. That it has also matured is indicated by the breadth 
of the subjects being studied with lattice methods. These include the traditional — QCD 
spectrum and matrix elements; the more statistical mechanical — rigorous theorems on finite 
size scaling; the more abstract — random surfaces; and the exotic — finite temperature baryon 
number violation. More concrete evidence for maturation is that lattice results are becoming 
useful to the rest of the particle physics community. For example, in his summary of B- 
physics, David Cassel noted that the bounds on the elements of the CKM-matrix which 
follow from B — B mixing depend on the values of and Bb, and that lattice results are 
now used as one of the estimates of these numbers. But perhaps the most important piece of 
evidence is that lattice studies are beginning to produce results with no unknown systematic 
errors. The best example is the calculation of the full QCD running coupling constant 
by the Fermilab group As I discuss below, this number can be directly compared to 
experiment. 

In summary, I would say that lattice studies are beginning to make themselves useful. 
The "beginning" in this claim is important — there is a very long way to go before we can 
calculate, say, the K ^ mr amplitudes from first principles. Thus this talk does not consist 
of results for a long list of matrix elements. Instead, I discuss a few topics in some detail. 
Because of lack of time, I concentrate entirely on results from lattice QCD (LQCD) which 
are relevant to phenomenology. 

When thinking about the progress that has been made, it is useful to keep in mind the 
questions that one ultimately wishes to answer using LQCD: 

1. Does QCD give rise to the observed spectrum of hadrons? Do these particles have the 
observed properties (charge radii, decay widths, etc)? 

2. Does the same theory also describe the perturbative, high-energy jet physics? 

3. What are the values of the hadronic matrix elements {fs, Bk, etc.) which are needed 
to determine the poorly known elements of the CKM matrix? 

4. What happens to QCD at finite temperature, and at finite density? 

5. What are the properties of "exotic" states in the spectrum, e.g. glueballs? 

6. What is the physics behind confinement and chiral symmetry breaking? 

The last question is the hardest, and, while it is important to keep thinking about it, there 
have not been any breakthroughs. The other questions are being addressed using numerical 
simulations, complemented by a variety of analytic calculations. Most recent progress has 
concerned the first four questions, and I will discuss only the first three. For additional details 
see the reviews of Mackenzie (heavy quarks), Petersson (finite temperature), Sachrajda (weak 
matrix elements) and Ukawa (spectrum) at LATTICE '92 0. 
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To set the stage I begin with a brief summary of LQCD, eschewing as many details as 
possible. The action of QCD is the sum of a gauge term, 

6 /■ 1 

^9 = ~^ T7^XlTr(G'^z/G'/,,,) , (1) 
g Jx iz 

with Gfj,i, the gluon field strength, and a quark part 

Jx q 

These expressions are valid in Euclidean space, where numerical lattice calculations are 
almost always done. QCD is put on a hypercubic lattice by placing the quarks on the sites, 
and gauge fields on the links which join these sites, in such a way that gauge invariance 
is maintained. The derivative in can be discretized in many ways, leading to different 
types of lattice fermions. The most common choices are Wilson and staggered fermions. The 
coupling constant g has been absorbed in the gauge fields, and appears only as an overall 
factor in Sg. It is conventional (see Eq. |T]) to use the combination P = 6/g'^ to specify g. In 
present simulations /? ~ 6, so that g^ ^ 1. The lattice spacing a is determined implicitly by 
the choice of g"^, as discussed below. 

The prototypical quantity of interest is the two-point correlator, e.g. 

C{t) = (j2[ho7,dit,x)] [J7o75&(0)]^ 

= sign(t) /^mfie-™^!*! (l + 0(e-(™s'-'"^)l*l)) . (2) 

At large Euclidean time, t, one picks out the lightest state (here the B meson). The contri- 
bution of excited states (beginning with the B') is suppressed exponentially. Thus one can 
just read off the mass, m^, while the amplitude of the exponential gives the decay constant 
(/b) up to kinematical factors. 

The expectation value in Eq. ^ indicates a functional integral over quarks, antiquarks 
and gluons weighted by exp{—Sg — Sq). Doing the quark and antiquark integrals, one obtains 

Cit) = - 1 [dA] e-^« [n,det(^ + m,)] ^ Tri^ol^Gdit, x; 0)7075^^(0; t, x)] . (3) 

Here [dA] is shorthand for the functional integral over lattice gauge fields, and the G are 
quark propagators Gq = [Ip + mq)~^ . The integral is normalized to give unity in the absence 
of the trace term. I have shown Eq. ^ diagrammatically in Fig. ||, where the lines are quark 
propagators, and the primed expectation value includes the determinant. 

To make the calculations numerically tractable, they must be done in a finite volume 
of A'^^ X Nt points. The functional integral is then of large but finite dimension, and can 
be done by Monte-Carlo methods. Similarly, the propagators are obtained by inverting a 
finite matrix. Available computer speed and memory limits the number of sites in this four 
dimensional world. In fact, it is speed that limits present simulations: the bottleneck is 
the inclusion of det(^ + rriq) in the measure. To make progress, many calculations use the 



2 




Figure 1: Diagrammatic representation of Eq. ^ 



"quenched" approximation in which the determinant in Eq. ^ is set equal to unity. This 
is equivalent to dropping all internal quark loops, so that there are only valence quarks. I 
discuss the consequences of this drastic approximation below. 

It is important to realize that numerical LQCD calculations which do not use the quenched 
approximation become exact when the lattice spacing vanishes, the volume goes to infinity, 
and the statistical errors vanish. They are thus non-perturbative calculations from first 
principles, which have errors that can be systematically reduced. They are not "model" cal- 
culations. On the other hand, calculations in the quenched approximation are more similar 
to those in a model. Certain physics is being left out, and one does not know a priori how 
large an error is being made. As I will show, present calculations suggest that the error is 
small, at least for a range of quantities. 

In the following three sections I discuss the spectrum, the calculation of and that of 
Jb- I will focus mainly on quenched results since these are more extensive. I will only make 
some brief comments on results for "full QCD" , i.e. QCD including the determinant in the 
measure. 

2 SPECTRUM 

I begin by discussing the status of calculations of the spectrum of light hadrons. I concentrate 
on m^r, rrip and {N refers to the nucleon), since these have the smallest errors. To give 
an idea of the rate of progress, I will compare with results from March 1990 0. The two 
and a half year interval since then is long enough to allow substantial changes. For example, 
computer power has increased as roughly CPU oc e^^^"^, so that, for a four dimensional theory 
such as QCD, the linear dimensions could have increased by a factor of ~ e^'^^^ = 1.8 in this 
period. In fact, the extra time has only partly been used in this way. 

To display the data I follow the "APE" group and plot m^r/mp versus (m^/mp)^ (see 
Figs. To understand the significance of these plots, recall the following. In a lattice 
calculation, we can dial the values of the quark masses. Ignoring for the moment the strange 
quark, and assuming degenerate up and down quarks, we then have a single light quark 
mass, nil, at our disposal. Each value of mi corresponds to a possible theory, each with 
different values for dimensionless mass ratios such as ml/rri^p, mp/rriN, fn/^p, m/^/m^, etc. 
We would like to fix mi using one of these ratios, and then predict the others. In practice, 
it is technically very difficult to do a simulation with small enough m;, and so we must 
extrapolate. The APE plot is one way of displaying how well this extrapolation agrees with 
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Figure 2: APE plots for quenched Wilson fermions: (a) March 1990 (b) November 1992 



the experimental masses. As mi varies from — >■ oo, ml/m^p varies monotonically from 
— >■ 1. Thus the theory maps out a curve in this plot, which we know must pass through 
the infinite mass point (m.^ = nip = 2m7v/3; shown by a square in the plots). The issue is 
whether this curve passes through the experimental point, indicated by a "?" on the plots. 

The solid lines if Figs. ^ are the predictions of two phenomenological models. The curve 
for light quark masses uses chiral perturbation theory, with certain assumptions, and is 
constrained to go through the physical point. The curve at heavier masses is based on a 
potential model. For more discussion see Ref. [^. 

What do we expect for the lattice results? There will be corrections to physical quantities 
which vanish in the continuum limit as powers of a (up to logarithms, which I will ignore), 
e.g. 

(mAr/mp)iatt = {mN/mp)cont[l + aA + 0{a^)] . (4) 

Thus, for finite a, the lattice curve should not pass through the experimental point. Similarly, 
if the physical size of the lattice [L = Nga) becomes too small the masses will be shifted 
from their infinite volume values. Thus, in addition to extrapolating in the quark mass, one 
must attempt to extrapolate both a — > and L —>■ (yo. 

What has happened in the last two years is that these extrapolations have become more 
reliable. I illustrate this with results using Wilson fermions in the quenched approximation 
(details of the data set are given in Table 0). Figure ^ shows the state-of-the-art in March 
1990. The upper two sets of points are from a large lattice spacing (a ~ 1/6 fm) with two 
different lattice sizes (L ^ 2 and 4fm). The results agree within errors, and I concluded that 
we knew the infinite volume curve for a = l/6fm with a precision of about 5%. This curve 
appeared not to pass through the physical point. The lower two sets of points are from a 
smaller lattice spacing (a ~ 1/10 fm). They suggested a downward shift in the data with 
decreasing lattice spacing, but it was difficult to draw definite conclusions given the large 
errors. 

The present results are shown in Fig. ^3. To avoid clutter, I show only data with 
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Table 1 : Parameters of lattices used to produce data shown in Fig. 
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a < l/9fm.lll These are all consistent with a single curve lying below that at a = 1/6 fm 
and close to the phenomenological predictions. The improvements in the results have come 
from using more lattices to approximate the functional integral, which reduces the statistical 
errors (see Table |lD. Furthermore, the decrease in a has been compensated by an increase 
in Ns, so that (with the exception of the APE results at /? = 6.3) the physical extent of 
the lattices exceeds L = 2fm. The results at /? = 5.7 imply that this is large enough to get 
within a few percent of the infinite volume results. 

The GFll group have used their results at a ~ 1/6, 1/9, 1/12 fm to do an extrapolation 
both to physical quark masses and to a = They find that, for a range of quantities, 
the results are consistent with their experimental values within ~ 5% errors. For example, 
rriN/mp = 1.284 (cf 1.222 expt.) and mA/mp = 1.627+J;[|^2 (cf 1.604 expt.). It should 
be recalled that a few years ago the quenched results for rriN/rnp were thought to be larger 
than the experimental value, while ttt-a — Tnj^ was smaller. What we have learned is that 
both the errors introduced by working in finite volume and at finite lattice spacing shift the 
curve in the APE plot upwards, so that one can easily be misled by results on small lattices. 

We seem, then, to know the spectrum of light hadrons in the quenched approximation 
fairly well, and it looks a lot like the observed spectrum. There is not, however, complete 
agreement on the numerical results. The QCDPAX collaboration, working at lattice spacings 
comparable to those in Fig. |^, finds significantly larger values for mi\i/mp at smaller m'^/m'^ 
PI . They suggest that the disagreement may be due to a contamination from excited states 
in the correlators (see Eq. |^). This disagreement will get cleared up in the next year or 
so. What is needed is operators which couple more strongly to the ground state, and less 
to excited states. There has been considerable improvement in such operators, for systems 
involving a heavy quark [^], but less so for light quark hadrons. 

Let me assume that the results of Fig. ^ are correct, so that the quenched spectrum 
does agree with the experiment to within ~ 5%. Does this imply that we can trust quenched 
calculations of other quantities to this accuracy? I do not think so. A priori, I would not 
have expected the quenched approximation to work so well, because it leaves out so much 
important physics. For example, the quenched rho cannot decay into two pions, unlike the 

^The new results at a = 1/6 fm ||^ confirm, with reduced errors, those in Fig. 
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physical rho, which might lead to a 10% underestimate of its mass |T^. Also, the pion cloud 
around the light hadrons is much different in the quenched approximation than in full QCD, 
which should shift rrip and mjv in different ways. While it is possible that these effects largely 
cancel for the spectrum, I see no reason for them to do so in other quantities. This argument 



can be made more concrete for the charge radii |TT[. Futhermore, as I mention below, there 



are reasons to think that the approach to the chiral limit in the quenched approximation is 
singular. 

My final remark on the quenched spectrum concerns the possibility of improving the 
approach to the continuum limit [|T2|. The gauge action is accurate to O(a^), but the Wilson 
fermion action has 0{a) corrections, as in Eq. It is possible to systematically "improve" 
the fermion action so that the corrections are successively reduced by powers of g"^ /Att. The 
idea is that by using a more complicated action one can work at a larger lattice spacing. The 
first results of this program are encouraging. Ref. finds that an improved action shifts 
the results at a = 1/6 fm downwards so as to agree with those at a ~ 1/10 fm, and that both 
agree with the "unimproved" results at a < 1/9 fm shown in Fig. |^. Ref. finds that 
the improved action makes no difference at a smaller lattice spacing, a ~ 1/14 fm. All this 
suggests that it may be possible use lattice spacings as large as 1/6 fm with an improved 
Wilson fermion action. Further evidence for this comes from lattice studies of charmonium 



T5| , p!6| . With staggered fermions, on the other hand, the prospects are less rosy. Although 
the corrections are of O(a^), and thus parametrically smaller than with Wilson fermions, 
they are in fact large at a = 1/6 fm |jl^ . 

Ultimately, we must repeat the calculation using full QCD. Much of the increase in 
computer time in the last few years has gone into simulations which include quark loops. 
Nevertheless, it is too early to discuss physical predictions, since it is not yet possible to 
reliably extrapolate ruq — >■ or a — 0. The limit L — cx) is, however, well understood and 
some interesting results have been obtained. The Kyoto-Tsukuba group finds good fits to the 
form m = moo(l + c/L^), and argue that this can be understood as due to the "squeezing" 
of the hadrons in the finite box [|l^.i The MILC collaboration finds that at L = 2.7 fm the 
finite volume effects are smaller than 2% ^^j. Both these results are consistent with the 



finite volume effects observed in the quenched approximation 



3 QMS FROM THE LATTICE 

It is straightforward, in principle, to extract from lattice calculations: 

1. Pick a value of (3 = G/g"^. 

2. Calculate a physically measurable quantity, e.g. rup or f.„. 

3. Compare lattice and physical values and extract a using, e.g.i 

("^p)latt = ("^p)phys X a (1 + 0(a)) . (5) 

The 0(a) terms are not included when extracting a. 

^This is not the asymptotic form: for large enough L the power law becomes an exponential ]l9t . 
§The values for a quoted above were obtained in this way, using a variety of physical quantities. 
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4. Convert from the lattice bare coupling constant g'^{a) to aj^(g = 7r/a) using pertur- 
bation theory. The result can then be evolved to other scales, e.g. Mz, using the 
renormalization group equation. 

5. Repeat for a variety of values of (3, and extrapolate to the continuum limit a = 0. In 
this way the 0(a) corrections in Eq. ^ are removed. 

If this program were to be carried through, then the lattice result for ajj^ would allow 
absolute predictions of jet cross sections, R{e~^e~), etc. (modulo the effects of hadronization). 
If these predictions were successful, it would demonstrate that QCD could simultaneously 
explain widely disparate phenomena occurring over a large range of mass scales. While such 
success has not yet been achieved, there has been considerable progress in the last year or 
so. I will explain the two major problems, and the extent to which they have been resolved. 
I then discuss the present results. 



3.1 Reliability of Perturbation Theory 

The fourth step in the program requires that perturbation theory be valid at the scale of 
the lattice cut-off, which is roughly ir/a in momentum space. On present lattices this ranges 
from 5 — 12 GeV (the values are given in Table |l]). It turns out that these values are not 
large enough for perturbation theory in the bare coupling constant to be accurate, because 
there are large higher order corrections. These are exemplified by the relation needed in step 



4 above 



^ = ^— -3.880 + 0(a) ; {a = g^Air) . (6) 

"ms(^/") ttiatt(a) 

Since aiatt ~ 1/13, the first order correction is large, and higher order terms are likely to be 
important. 

This problem has been understood by Lepage and Mackenzie |2^. The large corrections 



arise from fluctuations of the lattice gauge fields, and in particular from "tadpole" diagrams 
which are present on the lattice but absent in the continuum. The solution is to express 
perturbative results in terms of a "continuum- like" coupling constant, e.g. omom or a^^, 
with the scale q ~ n/a. When this is done the higher order coefficients are considerably 
reduced. This is similar to what happens in the continuum when one shifts from the MS to 
the MS scheme. 

Having re-expressed all perturbative expressions in terms of, say, a^jg, there remains the 
problem of finding the value of this coupling in terms of aiatt, since Eq. ^ is not reliable. 
One has to use a non-perturbative definition of coupling constant which automatically sums 



up the large contributions, and which is related to a^^ in a reliable way. One choice is p2 



31n(Trf/p) 1 1 n r ^/ n /^x 
ap = \ '-■ rT-\= 0.5 + O(«), 7 

where Up is the product of gauge links around an elementary square. One first determines 
ap from the numerical value of TiUp, and then converts this to Oms, which is then used in 
perturbative expressions. Lepage and Mackenzie find that the resulting numerical predictions 
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of lattice perturbation theory work well for quantities (such as small Wilson loops) which are 
dominated by short-distance perturbative contributions. This is true for lattice spacings as 
large as 1/9 fm, and perhaps 1/6 fm. Thus the determination of using Eq. |^ is reliable, 
with errors probably no larger than a few percent. 

3.2 Errors Introduced by the Quenched Approximation 

The dominant source of uncertainty in present calculations of is the use of the quenched 
approximation. The problem is the lack of a "physical" quenched theory to use in the 
comparison of step 3 above. This shows up in two ways. First, the value of a depends on 
the quantity chosen in the comparison: using mp gives one value, /tt another.B Second, the 
coupling constant that one obtains is for a theory with zero flavors, a^, and must somehow 
be related to the physical coupling constant a^jg. However, such a relationship involves 
non-perturbative physics, and so can only be determined by a calculation using full QCD! 
The best that we can hope for at present is a good estimate of the relationship between the 
couphngs. 

Recently, the FNAL group have made such an estimate for the coupling determined 
using the IP — IS* splitting in charmonium |l[]. The crucial simplifying feature is that 
charmonium is described reasonably well by a potential model, so one need only estimate 
the effects of quark loops on the potential itself. In outline, this is done as follows. Matching 
the lattice and continuum IP — IS" splittings makes the quenched and full QCD potentials 
similar at separations R ~ 0.5 fm. The potentials will, however, differ at smaller separations. 
We understand this difference at small enough R, where the Coulomb term dominates, i.e. 
V oc —ajj^{l/ R) I R. In the quenched approximation a varies more rapidly with R because of 
the absence of fermion loops, which means that the quenched potential is steeper. Assuming 
that this is true all the way out to 0.5 fm, where the potentials match, implies that the 
quenched potential must lie below that for full QCD at short distances. It follows that 
•^Ms ^ '^M^ short distances. A quantitative estimate is given in Ref. 

Unfortunately, there is no such simple way of estimating the effects of the quenched 
approximation on the values of extracted from the properties of light quark hadrons. 

3.3 Results 

Various physical quantities have been used to calculate and I collect the most accurate 
results in Table 0. I also give the "experimental" number obtained from an average of various 
perturbative QCD fits to data |2^. I quote the coupling at the scale M^, to allow comparison 



with Fig. 1 of the 1992 Review of Particle Properties (RPP) [^. The second row gives the 
results of the FNAL group, including the correction for quenching. The remaining results 
use the "string tension", o. This is the coefficient of the linear term in the heavy quark- 
antiquark potential: V{R) oR for R ^ oo. All groups use the "improved" perturbation 
theory explained above. I have not shown results from light hadron masses as they are less 
accurate. 



^Actually, since light hadron mass ratios are weU reproduced by the quenched approximation, the variation 
of a is small for such quantities. This need not be true in general. 
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Table 2: Results for in quenched and full QCD. Errors are statistical. 



Quantity Ref. a^^{Mz) ay^{Mz 



"Experiment'' 



m 0.1134(35) 
m 0.0790(9) 0.105(4) 



String tension Q 0.0796(3) 
String tension H 0.0801(9) 



The results for obtained using a have the smallest statistical errors. Indeed, it is a 

MS ° ' 

triumph of LQCD that the linear term in the quenched potential is very well established, for 
it is this term that causes confinement. The difficulty is that a is not a physical quantity. 
The "physical" value, ^/a = 0.44 GeV, is extracted from the potential needed to fit the cc 
and bb spectra. These systems do not, however, probe the region of the potential where 
the linear term dominates. Furthermore, at large R the full QCD potential fiattens out 
due to quark-pair creation, while the quenched potential continues its linear rise. Thus the 
extracted value of a is somewhat uncertain, and it is difficult to relate obtained using a 
to the full QCD coupling. Nevertheless, it seems to me possible to make a rough estimate, 
and it would certainly be interesting to try. 

The lattice prediction for ajj^ is slightly less than 2a below the RPP average. In his 
summary talk, Keith Ellis quoted an updated average, a^{Mz) = 0.119(4), which is almost 
3a higher than the lattice result. I think that this near agreement is a success of LQCD, 
and that it is too early to make anything out of the small discrepancy. For one thing, the 
dominant error in the lattice result is the uncertainty in the conversion from to ctj^, 
and this might be underestimated. It is important to realize, however, that this uncertainty 
will be gradually reduced as simulations of full QCD improve. 

I end this section with a general comment on the calculations. For the method to work, 
both non-perturbative and perturbative physics must be included in a single lattice simula- 
tion. At long distances the quarks are confined in hadrons, while at the scale of the lattice 
spacing their interactions must be described by perturbation theory. On a finite lattice, these 
two requirements pull in opposite directions. For example, if perturbation theory requires 
a < 1/9 fm, while finite volume effects require L > 2.7 fm (as may be true for light hadrons), 
the lattice must have Ng > 30, and so present lattices are barely large enough. This is 
another reason for using charmonium to determine a^g- The cc states are smaller than light 
hadrons, and it turns out that Ng = 24 is large enough to reduce the finite volume errors 
below 1% i. 

Nevertheless, it would be nice to extend the range of scales so as to provide a detailed 
test of the dependence on a. To accomplish this with present resources requires the use of 
multiple lattices having a range of sizes. Liischer et al. have proposed such a program, and 
done the calculation for SU(2) pure gauge theory p6| . 
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4 The B-meson decay constant 



One of the most important numbers that LQCD can provide phenomenologists is Jb- Anal- 
yses of the constraints due to K — K and B — B mixing find (for irtt > 140 GeV) two types of 
solutions for CP- violation in the CKM-matrix (see, e.g. Ref. [0). These are distinguished 
by whether is "small" (100 - 150 MeV) or "large" (160 - 340 MeV). In the former case, 
CP-violation in B ^ KsJ/ip is small; in the latter it is much larger. We would like the 
lattice to resolve this ambiguity. 

The calculation of fs is also interesting as a testing ground for the "heavy quark effective 
theory" (HQET) |^. The B-meson consists of heavy b quark, and a light quark {u or d). 
Imagine that we were free to vary the heavy quark mass, mg. The mass of the pseudoscalar 
meson (which I call mp to distiguish it from the physical B meson mass, m^), and its decay 
constant fp, would depend upon mg. As ttiq — > oo one can show that p9| 

<pp = fpy/rn^ 



a{mp) 


2//3o / 




= 0OO 1 



1+A+4+...V (8) 

nip nip J 

Here f3o is the first term in the QCD /5-function, and A and B are constants except for a 
weak logarithmic dependence on mp. The issue for HQET is the size of the 1/nip corrections 
at the B and D masses, for these might be indicative of the size of the corrections in other 
applications of HQET, e.g. B ^ D transitions. 

It is important to realize that, while the mq — >■ cxo limit simplifies the kinematics of the 
heavy quark, it does not simplify the dynamics of the light quark. In particular, the quenched 
approximation is no better in the heavy quark limit, nor is the dynamics more perturbative. 
Thus one needs a LQCD calculation here just as much as for light quark hadrons. 

What we would like to do is map out the curve of 0p versus 1/mp, and read off the 
values of fs and fn. Examples of present results are shown in Figs. |^, and in Table ^ I 
give the present numerical values for fs and fc. I also quote /p(static) = 0oo/ y/ms, which 
is the value of /p ignoring the 1/mp corrections in Eq. |^. Unfortunately, the situation is 
less straightforward than the figures imply, and I will spend the remainder of this section 
explaining and evaluating these results. 

There are three major causes of uncertainty in fp- The first concerns the overall scale. 
To extract a result in physical units we need to know the lattice spacing. As discussed above, 
the value we obtain depends on the physical quantity that we use, particularly at finite lattice 
spacing, fp is more sensitive to the uncertainty in a than are light hadron masses, since 
one is calculating a^^'^cpp rather than am. To illustrate this sensitivity, I include in Table 
^ results from Refs. [^, using two different determination of a. For the other results, 
the uncertainty is about 15%. Ultimately, to remove this uncertainty one must repeat the 
calculation with full QCD. 

The second problem concerns the isolation of the lightest state. In principle, calculating 
0P is straightforward. One simply studies the long time behavior of the two point function 
of the local axial current, Eq. 0, and reads off fpy/mp. In practice, to obtain a signal it is 
necessary to use extended operators, which couple more strongly to the lightest state than 
the local current. This is particularly important at large mg. With a less than optimal 
operator there are likely to be systematic errors introduced by contamination from excited 
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Figure 3: Results for (pp = rup^fmp: (a) from Ref. [^], (b) from Ref. |^ 



Table 3: Results for decays constants in the quenched approximation. The normahzation is such that 
/.„- = 132 MeV. Only statistical errors are shown; systematic errors are discussed in the text. 
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6.3 
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M{IP)-M{IS) 
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a 


265(10) 








5.74-6.26 


a 


230(22) 












256(78) 







states. Indeed, there are disagreements between results using different operators. This is 
illustrated in Table |^ by the variation in (static). It seems to me, however, that the 
operators of Ref. are close to optimal, and that the discrepancies will go away as other 
groups optimize their operators. 

The third and most difficult problem concerns putting very heavy quarks on the lattice. 
As the quark's mass increases, the ratio of its Compton wavelength to the lattice spacing, 
l/{mQa), decreases. For mqa > 1, its propagation through the lattice will be severely 
affected by lattice artifacts. There are, however, no such difficulties for an infinitely massive 
quark, for such a quark remains at rest both in the continuum and on the lattice |2^. Thus 
it seems that we are forced to interpolate between rriQ ~ 1/a and rriQ = oo. 

To illustrate the problems that this introduces, consider the situation about two years 
ago. The smallest lattice spacing was a = 1/10 fm, so that rrica ~ 0.75 and rnj, ^ 2.5. Typical 
results are those represented by squares in Fig. |^. The quark mass has been restricted to 
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satisfy mqa < 0.7 to avoid large artifacts (an arbitrary but reasonable choice for an upper 
bound), and thus lie to the right of the D meson line. Clearly it is difficult to convincingly 
interpolate using Eq. H: the variation in 0p is so large that one cannot truncate the 1/mp 
expansion. Nevertheless, if one assumes that the curvature is not too large, the fact that we 
know 0OO does allow a rough estimate of fs- 

There has been considerable progress in the last two years. The lattice spacing has been 
reduced to 1/17 fm, which allows one to use heavier quarks while keeping mqa < 1. This is 
illustrated by the remaining points in Fig. all of which have rriQa < 0.7. These points 
can now be fit in a more reliable way to the asymptotic form of Eq. ^. The curve shows 
such a fit, the results from which are given in Table ^ Other groups, however, find results in 
apparent disagreement. For example, the "uncorrected" points in Fig. ^ should agree with 
those in but instead are lower. This disagreement is, I suspect, mainly due to inadequate 
isolation of the lightest B meson by one or both groups. 



Another development has been the use an improved fermion action |Q. Since this has 
smaller 0{a) corrections, one should be able to work at larger values of mqa. I give the 
results in Table ^. 

Despite these improvements, it would be much better if we could work at any value of 
mqa, and simply map out the entire curve. There is a dispute about whether this is possible, 
and I will attempt to give a summary of the arguments. 

I begin by noting that the errors in 0p do not keep growing as mqa increases. This is 
because, if niq ^ Aqcd, the quark is non-relativistic. Its dynamics can then be expanded 
in powers of l/mq p5[ 



C = ip 



t 



^, (9) 



where ijj is a. two component field, and c is a perturbative coefficient. In the continuum, 
mi = m2 = mq. The discrete rotational symmetries are sufficient to ensure that a heavy 
lattice quark will be described by the same Lagrangian, except that neither mi nor m2 are 
equal to the bare lattice quark mass mq. Nevertheless, although the lattice heavy quarks 
have the wrong dynamics, this should only introduce errors of 0(Aqcd/'^q)- Thus, if a is 
small enough that mqa < 1 when the quark becomes non-relativistic, the errors will be small 
for all mqa. 

Strictly speaking, to bring the lattice Lagrangian into the form of Eq. one must 
perform a wavefuction renormalization on the heavy field, which changes the normalization 
of the bu axial current. This has been calculated in perturbation theory keeping only the 
"tadpole" diagrams which are thought to be dominant |]5B[. The result is illustrated in Fig. 

upon renormalization the "uncorrected" points are shifted upwards into those labeled 
"corrected" . With this modification, the curve for finite mqa is guaranteed to pass through 
the 0OO- Without it, the curve will bend over and eventually pass through the origin. In 
fact, the uncorrected points in Fig. ||b do appear to be fiattening out at the smallest values 
of 1/mp, although those of Fig. |^a do not. 

A second correction has also been applied in Fig. Pd. It is possible to partly account for 
the difference between mi and mq by an appropriate shift to smaller l/mq. The size of this 



shift has been calculated keeping "tadpole" diagrams [36]. This correction ensures that the 
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data approach (poo with a slope which is hnear in 1/mp. This slope will not, however, be 
correct since the 1 / m2 term in Eq. ^ has the wrong normalization. Nevertheless, in contrast 
to the "uncorrected" points, the shifted points fit well to the form of Eq. |^. The fit is shown 
by the curve, and the resulting values for decay constants are given in Table ^. 

The controversial issue is how to further modify the calculation so as to obtain the correct 
l/rrip terms. Kronfeld, Lepage and Mackenzie have suggested a program for removing the 



error [Q. The idea is to add new terms to the lattice action, and to choose the parameters 
so that one obtains the correct Lagrangian at ©(l/mg). This is essentially a complicated 
way of putting non-relativistic QCD on the lattice |^5|]. The difficult question in both cases 
is how to fix the parameters. 

Kronfeld et al propose that this can be done using perturbation theory. Maiani et al. 
argue, however, that non-perturbative contributions may be important Normally, such 
contributions are suppressed by powers of a, but here they are enhanced by factors of 1/a 



due to linear divergences. (A clear explanation of this effect is given in Ref. [^.) The 



result is that, with parameters fixed using perturbation theory, there will be errors in 0p 
which are of 0(Aqcd/'"^q), i-e. of the same size as the terms one is trying to calculate. The 
point of disagreement is whether these non-perturbative terms are large or small. If they 
are large, the only solution would be to fix parameters using non-perturbative normalization 
conditions. In general this would reduce predictive power. 

Clearly more work is needed to resolve this dispute. One way to do this is to push the 
tests of perturbation theory ||22| to the level at which non-perturbative terms show up. 

It is fortunate, however, that this uncertainty has only a small effect on fp and fo- This 
is shown by the good agreement in Table ^ despite the different methods being used. The 
results favor the "large" solution for fp- The only way this could change is if the systematic 
error due to the quenched approximation turns out to be large. 

We can also estimate the size of the 0(1/M) corrections to the heavy quark limit. Taking 
the data of Ref. |Q as an example, they are ~ 15% for /s, and ~ 45% for fo. These numbers 
increase further if one uses the larger values of 0oo found by Ref. 



5 OUTLOOK 

There are many interesting developments that I have not have had time or space to cover. 
One which impinges on much of the work discussed above concerns the accuracy of the 
quenched approximation In this approximation there are nine pseudo-Goldstone bosons, 
rather than the eight of QCD, the extra one being the rj' . This means that quenched particles 
have an rj' cloud, not present in full QCD. It turns out that this gives rise to singularities in 
the chiral limit, much as pion loops give singularities in the charge radii of pions and nucleons 



| TT| . For example, the quark condensate diverges. The implications of these divergences are 
not yet clear. The most optimistic view is that the effects of rj' loops are small as long as we 
work above a certain quark mass. This is supported by the absence of numerical evidence 
to date for the divergences. 

Unfortunately, the quenched approximation is likely to be with us for a number of years. 
At present, the fastest computers simulating LQCD are running at 1 — 10 GFlops. Even a 
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TeraFlop machine, such as that proposed by the TeraFlop collaboration [0, will focus on 
quenched lattices. These will have A^^ ~ 100, and should give definitive quenched results for 
a reasonable number of interesting quantities. The major problem with simulations of full 
QCD is that the CPU time scales as m~^°'^ with present algorithms. Clearly, it is crucial 
that effort go into improving these algorithms. 
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